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Abstract: Acceptance sampling is one of the essential areas of quality control. In a conventional 
environment, probability theory is used to study acceptance sampling plans. In some situations, 
it is not possible to apply conventional techniques due to vagueness in the values emerging from 
the complexities of processor measurement methods. There are two types of acceptance sampling 
plans: attribute and variable. One of the important elements in attribute acceptance sampling is 
the proportion of defective items. In some situations, this proportion is not a precise value, but 
vague. In this case, it is suitable to apply flexible techniques to study the fuzzy proportion. Fuzzy 
set theory is used to investigate such concepts. It is observed there is no research available to apply 
Birnbaum-Saunders distribution in fuzzy acceptance sampling. In this article, it is assumed that the 
proportion of defective items is fuzzy and follows the Birnbaum-Saunders distribution. A single 
acceptance sampling plan, based on binomial distribution, is used to design the fuzzy operating 
characteristic (FOC) curve. Results are illustrated with examples. One real-life example is also 
presented in the article. The results show the behavior of curves with different combinations of 
parameters of Birnbaum-Saunders distribution. The novelty of this study is to use the probability 
distribution function of Birnbaum-Saunders distribution as a proportion of defective items and find 
the acceptance probability in a fuzzy environment. This is an application of Birnbaum-Saunders 
distribution in fuzzy acceptance sampling. 


Keywords: fuzzy operating characteristic curve; fuzzy OC band; Birnbaum-Sunders distribution; 
single acceptance sampling plan 


1. Introduction 


An acceptance sampling plan is used to determine how many units can be selected from a lot, or 
consignment, and how many defective units are allowed in that sample. If the number of defective 
units is above the preset number of defective items, the lot is excluded. According to the rule of 
acceptance sampling, quality can be monitored by checking a few units from the whole lot. The 
plan that mentions guidelines for sampling and the associated criteria for accepting or rejecting a 
lot is called the acceptance sampling plan. This acceptance sampling plan can be implemented to 
check raw material, the material in a process or finished goods. An acceptance sampling plan can 
be classified as an attribute acceptance sampling plan and a variable acceptance sampling plan. An 
acceptance sampling plan can be classified with further attributes as a single sampling plan, double 
sampling plan, multiple sampling plans, and sequential sampling plan. An elementary acceptance 
sampling plan is a single sampling plan. In a single sampling plan, we select (n) units from the 
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entire lot. This consists of (N) units. After selection of n units they are examined; if the number of 
damaged units (d) is more than the specified number of defective items (c), the lot will be disallowed. 
Otherwise, it will be passed. The performance of any acceptance sampling plan can be judged by its 
operating characteristic (OC) curve. It determines how well an acceptance sampling plan distinguishes 
between good and bad lots. This OC curve has two parameters (n, c), where n is sample size and c 
is acceptance number. In an acceptance sampling plan, two groups are involved: the supplier and 
buyer. The supplier desires to avoid rejection of a good lot (producer’s risk) and the buyer tries to 
avert acceptance of a bad lot (consumer’s risk). In case a bad lot is accepted, it is the responsibility 
of the consumer. The producer’s risk is denoted by «. This is the probability of rejection of the lot 
having an average quality level (AQL). Similarly, the consumer’s risk is denoted with 8. This shows 
the probability of acceptance of the lot, having low quality (LQL) [1]. The proportion of defective 
items is denoted by p and treated as a precise number. However, in some situations, it is not possible 
to get the precise numerical value of p. Mostly this value is determined by the expert, based on his 
judgment. It is used to calculate fuzzy acceptance probability. Further, this fuzzy p value and fuzzy 
acceptance probability are used to design a fuzzy OC curve [1]. In the study presented in Reference [2], 
the authors suggested a double acceptance sampling plan, based on assumption that lifetime of the 
product follows a generalized logistic distribution with known shape parameters, and analyzed the 
operating characteristic curve to several ratios of the true median life to the specified life. In the study 
presented in Reference [3], the authors proposed the double sampling plan and specified the design 
parameters fulfilling both the producer’s and consumer’s risks at the same time for a stated reliability, 
in the form of the mean ratio to the specific life. Moreover, double sampling and group sampling plans 
are constructed using the two-point technique, with the assumption that the lifetime of the product 
follows the Birnbaum-Saunders distribution. In the study presented in Reference [4], the pioneer of 
fuzzy set theory gave scientific structure to study imprecise and ambiguous concepts that are based 
on human judgment; comprising verbal expressions, contentment degree and significance degree, 
that are often fuzzy. A linguistic variable consists of expressions in a natural language, but not the 
number. In Reference [5] the authors applied fuzzy set theory to help explain complex and not easy 
to express linguistic terms, in traditional measurable terms. In Reference [6] the authors proposed a 
single acceptance sampling plan with a fuzzy parameter and explained the single acceptance sampling 
plan with fuzzy probability theory. In Reference [7] the authors used the expression for the OC curve 
and various values to help accept or reject a lot for a particular number of defective items. Proficiency 
of different acceptance sampling plans can be assessed by using the OC curve. These OC curves 
are used to determine the producer’s risk, as well as the consumer’s risk [8]. In Reference [9] the 
authors suggested using acceptance sampling in the fuzzy environment using Poisson distribution. In 
Reference [10] the authors explored if N is large, then the defective items will have a fuzzy binomial 
distribution. In Reference [11] the authors applied parameters of the acceptance sampling plan, sample 
size n, and acceptance number c, in a fuzzy environment. Acceptance probabilities of two major 
discrete distributions were also derived. The multiple deferred sampling plans and characteristic 
curves were proposed—where (p) proportion of defective items was treated as a fuzzy number—and 
also proposed fuzzy OC curves with different combinations of parameters [12]. Multiple deferred 
acceptance sampling plans with inspection errors were proposed by the authors of Reference [13]. 
In the study presented in Reference [14], the authors investigated the inspection errors and their 
impact on a single acceptance sampling plan, when the proportion of defective items was not known 
exactly. In Reference [15] the authors proposed an acceptance sampling plan for geospatial data with 
uncertainty in the proportion of defective items. In Reference [16] the authors investigated a double 
acceptance sampling plan with the fuzzy parameter. Average outgoing quality (AOQ) and average total 
inspection (ATI) in a double acceptance sampling plan with the imprecise proportion of defective items 
were presented [17]. In Reference [18] the authors suggested the fuzzy parameter for quality interval 
acceptance sampling plan, applying Poisson distribution. The fuzzy double acceptance sampling 
plan for Poisson distribution was proposed by the authors in Reference [19]. In Reference [20] the 
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authors proposed an application of Weibull distribution in an acceptance sampling plan in the fuzzy 
environment and calculated fuzzy acceptance probabilities for different sample sizes using real-life 
data. In Reference [21] the authors proposed truncated life time, based on the Birnbaum-Saunders (BS) 
distribution. This distribution is used to define the number of stress cycles until failure of the material. 
In Reference [22] the authors applied the concept of the failure process of materials due to weariness, 
to design the BS distribution. Estimation of parameters based on crack length data was proposed in 
Reference [23]. In Reference [24] the authors presented a literature review of the BS distribution and 
discussed in detail the importance of this distribution and its application in different fields. In this 
study [25], they developed an acceptance sampling plan using the BS distribution to get the minimum 
sample size, n. 

The aim of this article was to apply a single acceptance sampling plan when data were fuzzy and 
the proportion of defective items followed the BS distribution. According to the best of our knowledge, 
there is no work on the fuzzy plan using the BS distribution in the literature. In this paper, we will 
develop the fuzzy sampling plan using this distribution. The application of the proposed sampling 
will be given with the aid of a real example. 


2. Materials and Methods 


Design of Proposed Plan 


Probability distribution function (Pdf) of BS distribution 


Lif 4 
Fraa) = 0(22(5)), 0< t<o,A>0 (1) 
where « is the shape parameter and A is the scale parameter, ®(.) is the standard normal cumulative 


function and ¢(t/A) = t — T . It can be shown that the median of the BS distribution is equal to 


the scale parameter and the mean of the BS distribution is 
w= A(1+0?/2) (2) 


Here we write the assumptions for the BS distribution. 
Let to = apo; a be called the termination ratio. The cumulative distribution function (Cdf) given 
in Equation (5) can be rewritten as 


X a? 
mana) =o(2e( 2/2) o 


aX u/ po 


The acceptance probability 

According to [26], the acceptance probability of sampling plans can be obtained by using the 
binomial distribution. The lot acceptance probability of a lot in a single acceptance sampling plan 
(SASP) case is given as 


n 


n . =2 

L(p) = pa Jap (4) 
i=0 

The proportion of defective items in the fuzzy form. 

According to the equation proposed by the authors of Reference [27]. 


pK = (K, bo + K,b3+ K, b4 + K).px = pK|«|,qx = gK|aæ], pK + 9x =1 (5) 


b; = 4; — a,i = 2, 3, 4 and K = [0,1 — bag]. 
x-cut of pK 
pK(«) = (K + (bz +K -— K)a, b3 + K+ (b3 = b4)a) (6) 
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a-cut of pK at x =0 
pK(0) = (K,b4 + K) 


where p is the Fr(to, &, À) in Equation (4). 
Fuzzy acceptance probability 
According to Reference [11], the fuzzy acceptance probability can be calculated as 





mat ( {eta 'lpe pageant osas: (7) 
Pr & = [Pry Prr] 
: n i 
pua=mind (Jeg ‘nepaacaa} (8) 
pra mast È : pc paacanh (9) 
The fuzzy acceptance probability when the number of defective items, c = 0, and « = 0 

ie = n ‘ae n 

Px (0)[0] = (1— Pr” [a]) , (1 -— Pr" la]) (10) 


DK’ =Kpx" =b4 + K 


The fuzzy acceptance probability based when the number of defective items, c = 1 and « = 0 
_ n = n—1 JS n pE n—1 
Px(1)[0] = (1—Pe“la]) +n (1—Pe“la]) (1 PK" [a]) +n (1 — Pral) (11) 


Here px" [a], px“ [a] are calculated using CDF of the BS distribution 
The design for a single acceptance sampling plan (ASP) to generate a fuzzy operating 
characteristic curve (FOC) 


Step 1. The sample size for a lot is n. 

Step 2. Specify the acceptance number (or action limit) c for a sample and the experiment time to. 

Step 3. Perform the experiment for the sample size n and record the number of failures for a sample. 

Step 4. Accept the lot if at most c failures are observed in the sample. Truncate the experiment and 
reject the lot if more than c failures are observed in the sample. 

Step 5. Calculation of fuzzy p (proportion of defective items) using Equation (6). 

Step 6. Calculation of fuzzy acceptance probability using Equations (10) and (11). 

Step 7. Design of fuzzy OC curve (FOC) includes k and the acceptance probability, where k is the 
transformation of the fuzzy proportion of defective items. 


The advantage of the fuzzy OC curve is that it is flexible and can be applied when the proportion 
of defective items is fuzzy. Secondly, the width of the fuzzy OC curve indicates the quality. Where 
the wider the width, the lesser the quality, and vice versa. In this study, the width of the fuzzy OC 
curve is influenced by the mean ratio. When the mean ratio is higher, the width of the band decreases. 
When the mean ratio is lower, the width increases. The advantage of this approach is that it can be 
applied to study any fuzzy data, which follows the BS distribution. This approach is more flexible 
than conventional p because it considers intermediate values of the fuzzy curve. 

The fuzzy proportion of defective item k at « = 0 is denoted by px|0] and the fuzzy acceptance 
probability as Px(0)[0]. Values for sample size n = 5 and acceptance number c = 0, will therefore be 
(0.00, 0.001), and (0.95, 0.96), respectively at k = 0.01. Similarly, values of proportion and acceptance 
probability for sample size n = 5 and acceptance number c = 0, will be (0.052, 0.054) and (0.77, 0.78), 
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respectively at k = 0.01. Furthermore, the fuzzy proportion of defective item k at « = 0 is denoted 
by px [0] and the fuzzy acceptance probability as Px(1)[0]. These values for sample size n = 5 and 
acceptance number c = 1, will be the proportion of defective items (0.01, 0.029), and the acceptance 
probability (0.99, 0.995) at K = 0.01. 


3. Real Life Example 


In this section, we will discuss the application of the proposed sampling plan using real data 
selected from [28] and [29]. As mentioned above, the BS distribution is also known as the fatigue life 
distribution. It is used extensively in reliability applications to model failure times. The BS distribution 
is used in circumstances where occurring of events is independent of each other, from one cycle 
to another cycle, with same random distribution [29]. In this study, the failure life data given in 
Reference [28] is used. The authors of Reference [28] found that the data follow the BS distribution. 
The fatigue life data of aluminum coupons having n = 101 observations are shown in Table 1. 


Table 1. Fatigue data of aluminum in hours. 


o|o |o | 1] 8) oe] ioe o 17] 18 | 18 | 8] 10 


162 | 163 | 163 | 164 | 166 | 166 | 168 | 170 | 174 | 16 faz 





We assume that data follows the BS distribution, the proportion of defective items p is fuzzy and 
the shape parameter is taken as the trapezoidal fuzzy number, & = (0.15, 0.16, 0.17, 0.18). When the 
actual mean is fp = 134, the termination ratio a = 0.5 is then truncated, time will be t = 67 for c = 0, 
and at oF = 1, the proportion of defective items is p = 0.0211673. The acceptance probabilities are 
calculated by using Equations (9) and (10), for different sample sizes n = (5,25,75,100) and K = (0.0, 
0.01,0.02, 0.03, 0.04, 0.05). The fuzzy OC curve is designed using fuzzy p values and fuzzy acceptance 
probabilities for c = 0. Similarly, when ig = 134, termination ratio a = 0.67 is then truncated, time will 
be t = 89.7 for c = 1. In this case, acceptance number c = 1, and ® = 1, p = 0.01923. The acceptance 
probabilities are calculated using Equation (11) for different sample sizes n = (5,25,75,100) and K = 
(0.0, 0.01, 0.02, 0.03, 0.04, 0.05). The fuzzy OC curve is developed using p values and acceptance 
probabilities for c = 1. Fuzzy acceptance probability for Birnbaum-Saunders distribution is presented 
in Tables 2 and 3 using real life data and their respective fuzzy OC curves are shown in Figures 1-3. 
Entire calculations and graphs were completed using R software and codes were given in Appendix A. 
Acceptance probability is influenced by mean ratio, when mean ratio increases it reduces Uuncertainty 
and bandwidth of fuzzy OC curve become narrow while decreasing mean ration increases the width 
of fuzzy OC curve. The fuzzy OC curves show more convexity when sample size n increases. Fuzzy 
OC curve with c = 0 shows less uncertainty than c = 1. We presented acceptance probabilities and 
fuzzy OC curves with c = 0, it is almost equal to conventional OC curve. The fuzzy OC curve for c = 0 
and c = 1 is more convex at large sample size as compared to small sample size. 


Table 2. Fatigue life data of aluminum coupons. 


Px(V[0], n=5,c=0, Pr(1)[0], n=25,c=0, PKIO] n=75,c=0, Px(1)[0], n = 100, 





K Pxl0] Hay Hey Hay c=0,h =1 
0.00 (0.00, 0.001] [0.98, 1.00] [0.99, 1.00] [0.99, 1.00] [0.996, 1.00] 
0.01 (0.01, 0.011] [0.95, 0.96] (0.81, 0.82] [0.73, 0.76] [0.604, 0.620] 
0.02 (0.02, 0.022] [0.93, 0.94] [0.65, 0.66] [0.53, 0.54] [0.3640, 0.372] 
0.03 (0.03, 0.041] [0.85, 0.86] [0.53, 0.54] [0.40, 0.43] [0.213, 0.219] 
0.04 (0.04, 0.051] (0.81, 0.83] (0.41, 0.43] (0.29, 0.31] [0.125, 0.129] 
0.05 [0.052, 0.054] [0.77, 0.78] [0.35, 0.36] (0.21, 0.23] [0.076, 0.077] 
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Table 3. The fuzzy acceptance probability of Birnbaum-Saunders distribution for c = 0. 


Px(D[0],n=5,c=1, Pr([0],n=25,c=1, Px(1)[0],n=75,c=1, Px()[0],n=100,c=1, 
H H H H 











K Bx 01 t - - a 
mal m=! m=! m=! 
0.00  [0.00, 0.019] [1.000, 1.00] [1.0000, 1.00] [1.00, 1.00] [0.99, 1.00] 
0.01  [0.01, 0.029] [0.99, 0.995] [0.979, 0.97] [0.827, 0.82] [0.82, 0.82] 
0.02  [0.02, 0.039] [0.94, 0.96] [0.99, 0.949] [0.55, 0.556] [0.65, 0.660] 
0.03 [0.03, 0.077] [0.93, 0.95] [0.82, 0.826] [0.338, 0.338] [0.63, 0.54] 
0.04  [0.04, 0.050] [0.91, 0.92] [0.73, 0.730] [0.190, 0.190] 0.44, 0.42] 
0.05 [0.05, 0.067] [0.89, 0.900] [0.720, 0.7202] [0.160, 0.160] [0.35, 0.37] 
(a) fuzzy OC curve (b) fuzzy OC curve 
O 
0 
O 
È o 
O 
T 
S 
0.00 0.02 0.04 0.00 0.02 0.04 
k k 
(c) fuzzy OC curve (d) fuzzy OC curve 
0 
O 
fy 
O 
Q 
S 
0.00 0.02 0.04 0.00 0.02 0.04 
k k 


Figure 1. The fuzzy operating characteristic (OC) Curve of the Birnbaum-Saunders distribution at c = 0 
(a) n = 20, (b) n = 40, (c) n = 70, and (d) n = 100. 


(a) fuzzyOC curve 


a=0.16,a=0.5 


) Qg=0-17,a=05 
5 x 0.18,a=0.5 
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g 2V,c=1 £ 
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044 oe ŒB 10 


OH 088 





Qo 


Figure 2. The fuzzy operating characteristic (OC) curve of Birnbaum-Saunders distribution at c = 1, 
mean ratio = 0.5, (a) n = 20, (b) n = 40, (c) n = 70, and (d) n = 100. 
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Figure 3. The fuzzy operating characteristic (OC) curve of Birnbaum-Saunders distribution at c = 1, 
mean ratio = 0.5, (a) n = 5, (b) n = 15, (c) n = 30, and (d) n = 50. 


4. Conclusions 


Acceptance sampling is one of the important aspects of statistical quality control. When the data 
follow the Birnbaum-Saunders distribution and the proportion of defective items is fuzzy, acceptance 
probability and the OC curve can be presented in a fuzzy form. In this article, the fuzzy OC curve of the 
Birnbaum-Saunders distribution is presented in a single acceptance sampling plan, using the binomial 
distribution. The fuzzy OC curve has a band with two bounds, lower and upper. The width of the 
band depends upon the uncertainty in the proportion of defective items in the fuzzy environment. Less 
uncertainty will give a narrow width. The fuzzy OC curves also show more convexity at a large sample 
size. The mean ratio in the Birnbaum-Saunders distribution is another important factor in quality. 
Here, a lower value of the mean ratio causes the width of the band of the fuzzy OC curve to increase. 
This indicates more uncertainty. The advantage of this approach is that it can be used to calculate the 
proportion of defective items when fuzzy data follows a Birnbaum-Saunders distribution, because 
mostly we assume the value of the proportion of defective items without using any distribution. 
Secondly, the fuzzy acceptance probability based on the Birnbaum-Saunders distribution is calculated. 
The fuzzy OC curve of the Birnbaum-Saunders distribution is constructed based on fuzzy p and the 
acceptance probability. The OC curve is more convex at large sample sizes, as compared to small 
sample sizes. It was concluded that when data followed the Birnbaum-Saunders distribution, this 
proposed approach was suitable to calculate the proportion (p), the acceptance probability, and the OC 
curve in both conventional and fuzzy form. In the future, we will apply the same concept to group 
acceptance sampling and chain acceptance sampling in a fuzzy environment. 
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Abbreviations 

OC operating characteristic 

OC curve operating characteristic curve 
(FOC) curve fuzzy operating characteristic curve 
AQL average quality level 

LOL low quality level 

pdf probability density function. 

cdf Cumulative distribution function 
SASP single acceptance sampling plan 
BS Birnbaum-Saunders 

Appendix A 

R codes 


#When n = 5, 20, 30, 30,c = 0 

rm(list=ls ()) 

windows () 

par(mfrow=c(2,2)) 

a= 0.5 #For c = 0 at t = 67 

alpha = c(0.15, 0.16, 0.17, 0.18) 

K = seq(0, 0.05, 0.01) 

x = a*((1 + alpha*2)/2) #Here b is teated as Alpha 

y = 0.5# 1, 2,3, 4, 5, 6 values of ratio of 

X = c((1/alpha)*(sqrt(x/y)—sqrt(y/x))) 

FX = pnorm(X, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE) 

FX 

p = FX 

p = 2.470246e-16 

K = seq(0,0.05,0.01) 

W=p+K 

p = 0.226627400 

W=p+K 

#p = 0.226627400 

p = 0.019 

W=p+kK 

K = seq(0,0.05,0.01) 

c = 0, a = 0.5) and (c = 1 when a = 0.67) for a = 0.5, p = 2.470246e-16, for a = 0.67 p = 0.01923 
B = dbinom(0,10,K) # B = (1—K)*5 When n = 5, c = 0 (1) 

A = dbinom(0,10, W)# A = (1—(K+p))*5#When n = 5, c = 0 (2) 

data.frame(A,B) 

data.frame(K,W,A,B) 

#B = (1-K)*5 # THIS WILL GIVE US UPPER BAND HIGHER PROBABILITY 

#A = (1-(K+p))5#THIS WILL GIVE US LOWER BAND LOWER PROBABILITY 
plot(K,A,type = “1”, col = “red”, xlab = “K”, ylab = “Pa “, main = “fuzzy OC curve”) 
par(new = TRUE) 

plot(W,B,type = “T”, col = “blue”, xlab = “k”, ylab = “ “, main = “”) 
legend(“topright”,c(expression(paste(alpha==0.15,”,”,a==0.67)),expression(paste(alpha==0.16,”,”,a== 
0.67)),expression(paste(alpha==0.17,”,”,a==0.1)),expression(paste(alpha==0.18,”,”,a==0.67)),expression 
(paste(n==5,”,’”",c==0)))) 
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